Leveraging opposition-based learning for solar photovoltaic model parameter estimation with exponential distribution optimization algorithm

Given the multi-model and nonlinear characteristics of photovoltaic (PV) models, parameter extraction presents a challenging problem. This challenge is exacerbated by the propensity of conventional algorithms to get trapped in local optima due to the complex nature of the problem. Accurate parameter estimation, nonetheless, is crucial due to its significant impact on the PV system’s performance, influencing both current and energy production. While traditional methods have provided reasonable results for PV model variables, they often require extensive computational resources, which impacts precision and robustness and results in many fitness evaluations. To address this problem, this paper presents an improved algorithm for PV parameter extraction, leveraging the opposition-based exponential distribution optimizer (OBEDO). The OBEDO method, equipped with opposition-based learning, provides an enhanced exploration capability and efficient exploitation of the search space, helping to mitigate the risk of entrapment in local optima. The proposed OBEDO algorithm is rigorously verified against state-of-the-art algorithms across various PV models, including single-diode, double-diode, three-diode, and photovoltaic module models. Practical and statistical results reveal that the OBEDO performs better than other algorithms in estimating parameters, demonstrating superior convergence speed, reliability, and accuracy. Moreover, the performance of the proposed algorithm is assessed using several case studies, further reinforcing its effectiveness. Therefore, the OBEDO, with its advantages in terms of computational efficiency and robustness, emerges as a promising solution for photovoltaic model parameter identification, making a significant contribution to enhancing the performance of PV systems.

Addressing climate change and shaping effective energy policies have become urgent global priorities.In this context, the value of photovoltaic (PV) power generation cannot be overstated.It offers a pathway to harness electricity from solar radiation, all while bypassing greenhouse gas emissions 1,2 .The usage of PV systems has witnessed a substantial rise, despite their relative costliness, as part of the solution to these global issues.Recent research has concentrated on crafting precise models to maximize the productivity of renewable energy (RE) sources.These endeavours aim to provide viable alternatives to traditional fossil fuels, whose usage is tied to environmental pollution.Among various renewable sources, solar energy stands out.Its capacity to generate electricity through PV power is a green alternative that causes no environmental harm.PV systems, pivotal in renewable energy development, transform solar energy into electrical power [3][4][5] .However, they face challenges due to environmental influences like dust, weather changes, and temperature variations.These elements can decrease the efficiency of solar cells, making the need for accurate PV models crucial.Such models can help optimize the • Incorporating the OBL technique is designed to enhance the quality of positions by utilizing a history of prior positions.• The comparison of OBEDO with other recognized algorithms using various PV models.
• Through comprehensive experimental results and statistical analyses, this study demonstrates the superior performance of OBEDO.
The rest of this document is organized as follows: Section "Photovoltaic modelling and problem formulation" delves into the various photovoltaic models, meticulously explaining their mathematical aspects, and explores the crafting of the objective function.Section "Proposed opposition-based exponential distribution optimizer" provides a succinct overview of the exponential distribution optimizer, thoroughly elaborating on the structure of the suggested OBEDO methodology.Section "Results and discussions" examines the results of various case studies, providing a profound analysis of the experiments.The final section wraps up the study, offering future research directions.

Photovoltaic modelling and problem formulation
This section details the various photovoltaic modes and their respective mathematical modelling.Furthermore, the objective function construction is also deliberated.

Photovoltaic modelling
Photovoltaic devices translate sunlight straight into electricity utilizing the photovoltaic effect.Accurately depicting these mechanisms is crucial to predicting their performance, designing systems, and undertaking thorough analysis.The three main models applied in PV system representation are SDM, DDM, and TDM 16,[81][82][83] .Typically, a cell is represented through a singular current source, denoted as I ph .The photocurrent, I ph , is reliant on the intensity of solar radiation.
The single-diode model is the most elementary and prevalently employed model for assessing PV systems.This model encompasses one semiconductor junction and the subsequent photovoltaic effect.Despite its simplicity, this model can effectively depict the function of a solar cell across diverse operational conditions.In the framework of the SDM, we engage a sole non-ideal diode, integrated into parallel alignment with the current source, a scenario visually represented in Fig. 1.As was broached earlier, R sh symbolizes the parallel resistance, R se stands for the series resistance, and the current traversing through the diode is designated as I d .However, a contrasting approach is adopted in the DDM, wherein a pair of imperfect diodes are collaboratively connected in a parallel manner to one unique current source, a setup graphically elucidated in Fig. 2. In this configuration, I d1 and I d2 denote the currents that flow through the first and second diodes, respectively, with I d2 in particular, corresponding to losses experienced within the space charge region.The fundamental mathematical formula for the single-diode model is as follows 84,85 .
The diode current I d is deduced using the Shockley diode formula, and the leakage current I sh is attributed to the shunt resistance, which reflects power losses in the PV cell.By including the expressions for I d and I sh , Eq. ( 1) is expanded as follows.
(1) I pv = I ph − I d − I sh ( The HSOA algorithm, which improved the basic SOA's optimization capability with three modified strategies, has been developed Table 1.Literature study. where I sd represents the diode saturation current, V pv denotes the PV cell output voltage, I pv denotes the PV cell output current, n denotes the ideality factor, V t is the thermal voltage, k denotes Boltzmann's constant, q denotes the electron charge, and T denotes the cell temperature.Examining Eq. ( 2) reveals five undetermined variables: I ph , n , I sd , R se , and R sh .Accurate estimation of these parameters from the I-V characteristic of the PV cell is crucial for successful PV modelling.However, the SDM does not consider losses in the depletion region.This has led to the proposition of the DDM by researchers who have asserted that the DDM yields more precise results.The DDM is a more intricate model accounting for recombination losses in both the depletion and quasi-neutral regions of the solar cell, phenomena not represented in the SDM.The equivalent circuit of the DDM is shown in Fig. 2. Essentially, it embodies two single-diode models functioning in parallel, hence "double-diode".The equation for the DDM becomes more sophisticated, incorporating two diode currents ( I d1 and I d2 ) is as follows 86 .
By replacing the expression for I d1 and I d2 using the Shockley equation, Eq. ( 4) is rewritten as follows.
where I sd1 and I sd2 are the saturated currents of diode 1 and diode 2, respectively, and n 1 and n 2 are the ideality factors of diode 1 and diode 2, respectively.
(3)  www.nature.com/scientificreports/Examining Eq. ( 5) reveals seven undetermined variables: I ph , n 1 , n 2 , I sd1 , I sd2 , R se , and R sh .Accurate estima- tion of these parameters from the I-V characteristic of the PV cell is crucial for successful PV modelling.Each diode current in DDM possesses its ideality factor, which enhances the precision of this model and its complexity compared to the single-diode model.The TDM incorporates a third diode to signify the recombination losses in the space-charge region, consequently achieving an even more precise representation of the actual behaviour of the solar cell.The equivalent circuit of the TDM is shown in Fig. 3.However, it also significantly elevates the complexity of the model, with the formula involving three diode currents (I d1 , I d2 and I d3 ) is as follows 87,88 .
By replacing the expression for I d1 , I d2 and I d3 using the Shockley equation, Eq. ( 6) is rewritten as follows.
where I sd1 , I sd2 and I sd3 are the saturated currents of diode 1, diode 2, and diode 3, respectively, and n 1 , n 2 , and n 3 are the ideality factors of diode 1, diode 2, and diode 3, respectively.Examining Eq. ( 7) reveals nine unde- termined variables: I ph , n 1 , n 2 , n 3 , I sd1 , I sd2 , I sd3 , R se , and R sh .Accurate estimation of these parameters from the I-V characteristic of the PV cell is crucial for successful PV modelling.Each diode current in the TDM also has its ideality factor.Despite being the most precise model, its complexity can lead to more challenging parameter extraction and increased computational burden.
These models support the comprehension and prediction of PV cells' behaviour under diverse illumination and temperature conditions.The SDM is ideal for applications where simplicity and computational speed are key, whereas the DDM and TDM offer greater accuracy and detail and are more suited for in-depth research and comprehensive analysis of PV cell operations.Building a PV module around the principles of the SDM involves integrating SDM-based PV cells arranged in series and parallel networks as key components.A PV module is structured using multiple parallel strings of PV cells ( N sh ), each containing an equivalent number of series- aligned PV cells ( N se ).The I-V characteristics of the module are derived using Eq. ( 8).

Problem formulation
The PV system, with its nonlinear, intuitive, and transcendental traits, makes the PV cells/modules an appealing choice for optimization tasks.An error function that assesses the precision of parameter estimates can feasibly be established by synthesizing Eqs.(2, 5, 7, and 8).When employing optimization algorithms for this task, creating an error or fitness function is a must; this function must be minimized to secure the best parameter estimation values.The selection of an error function is crucial since it has a substantial impact on the general effectiveness of the final model.Root mean square error (RMSE), as outlined in Eq. ( 12), is used as the objective function, providing a comprehensive overview of the performance of the estimated model across all characteristics.
The corresponding functions for the SDM, DDM, and TDM, along with the solution vector, denoted as x , for each PV model, are presented as follows.

Exponential distribution optimizer
The exponential distribution optimizer (EDO) is a metaheuristic algorithm rooted in mathematics that addresses intricate optimization problems 79 .The algorithm is founded on Exponential Distribution, a specific form of the continuous probability distribution that models the time between events in a Poisson point process.EDO effectively solves diverse optimization problems, including continuous, linear, nonlinear, and constrained problems.
The idea of EDO is to employ a population of individuals (solutions) to explore the search space.Each solution represents a point in the search space, and the objective function value determines the quality of the solution.
The algorithm exploits the balance between exploration (global search) and exploitation (local search) to direct the search procedure towards the global optima.In the phase of exploitation, the EDO leverages three fundamental constituents intrinsic to the exponential probability distribution (EPD).These constituents encompass the memoryless attribute, the directive solution, and the exponential dispersion, all essential to retain the contemporaneity of the innovative solution.Conversely, during the phase of exploration, a model optimized based on two solutions, both derived from the ED present within the initial population, is selected.Concurrently, the arithmetic average solution is utilized to effectuate an update to the promising solution.
The main inspiration of EDO, the ED, is a continuous probability distribution that describes the time between events in a Poisson point process, i.e., a process in which events occur continuously and independently at a constant average rate.The probability density function (PDF) of an exponential distribution is given by: where = 1 Mean .The ED has the property of being memoryless.In the context of the EDO algorithm, this implies that the quality of a solution does not depend on how the solution was obtained in previous iterations.In adherence to the memoryless property, there is no retention or consideration of the prior history of the solutions.This is due to the independence of past failures, rendering them empty of any impact on subsequent outcomes.Updated solutions are duplicated into the memoryless matrix to emulate the memoryless attribute intrinsic to the ED, irrespective of their fitness levels.This is premised on the understanding that historical data doesn't influence future developments.Consequently, the memoryless matrix becomes a repository for two categories of solutions: successful ones and those that fall short.Initially, the memoryless matrix is assigned a value identical to the original population, designated as X w .Like all metaheuristic algorithms, the initialization phase is done random population solution, i.e., the populations are randomly aligned, and the solutions are generated randomly using the random populations.Equation ( 14) can be employed to stochastically generate each exponential random parameter that is part of the candidate ED within the domain of the problem.
where lb and ub denote the lower and upper bounds of the control vectors, rand denotes the uniform random number varies between [0,1], N p signifies the population size, and dim denotes the problem dimension.
The exploitation stage leverages numerous features of the ED model, including the memoryless attribute, mean, standard deviation, and exponential rate.Furthermore, a directive solution is utilized to steer the exploration phase towards the global optima.The vicinity of a proficient solution often holds the potential for identifying the global optima.This is why several procedures delve into the search bounds near effective solutions by drawing (10)   www.nature.com/scientificreports/ in the lesser-performing ones.Hence, the quest for the global optima is centred on the guiding solution X g .The X g is characterized as the average of the top three solutions from an organized population, designed as follows: where t represents the current iterations, the guiding solution is chosen over the best solution because it incre- mentally leads the solutions towards the optimum one.Even if the best solution is ensnared in local optimum, all other outcomes persistently gravitate towards this best solution.The area surrounding an effective solution often harbours the potential for uncovering the global optimum.Consequently, numerous algorithms take advantage of the search space near high-performing solutions by drawing in the less successful ones.The exploitation stage of the EDO practices several features of the ED model, such as the memoryless feature, mean, and standard variance, to update the solution as follows.
where ML t i denotes the i th solution of the memoryless matrix, ϕ signifies the uniform random number between [0, 1], b and a denote adaptive variables and f denotes the random number in the range of [− 1, 1].The expression for exponential variance is provided as follows.
The exponential mean µ is determined as the mean value between the X t g and the i th memoryless parameter, which could either be a winning or a losing solution.
The algorithm's exploration stage pinpoints areas within the search bound that are considered likely to contain the globally optimal results.The optimization framework for the EDO's exploration stage is constructed around two successful solutions from the initial agents in adherence to the ED.Following this, the solution is updated using Eqs.(23-27).
where M t denotes the mean solutions obtained in the actual population, c denotes the tunable parameter signify- ing the ratio of information shared among Y 2 and Y 1 to the current solution, T denotes the maximum number of iterations, and X t w rand1 and X t w rand2 denote the winners concerning the randomly selected ED from the population.An approach based on a greedy method is utilized for the solution update obtained during the exploration and exploitation stages within the initial population.Any updated solution is integrated into the original population only if it meets the criteria of being considered good.The pseudocode of the EDO algorithm exists in Algorithm 1. (15) Initialize the parameters, such as , , bounds, etc. Initialize the random population and generate the random population solutions using Eq.14. Construct the memory matrix.

While ( < ) do
Define the matrix of size Compute using Eq. 15 Define all the adaptive parameters , , , , if ( < 0.5) Update the position using Eq.16 (Exploitation phase) else Update the position using Eq.23 (Exploration phase)

End if
Copy the matrix to the memoryless matrix.Update the solution comparing the current fitness and updated fitness.End while Return: Best position and the respective solution.
Algorithm 1. Pseudocode of the EDO algorithm.

Opposition-based exponential distribution optimizer
The mutation operations in the EDO algorithm are based on the ED parameters and two control parameters ( a and b ) derived from a uniformly distributed random variable f.This ensures diversity in the population and facilitates exploration and exploitation of the search space.Opposition-based learning (OBL) is an additional strategy that can be incorporated into EDO to enhance performance 80,89 .The foundational concept of OBL is based on De Morgan's laws in Boolean logic, which state that the complement of a conjunction is the disjunction of the complements, and the complement of a disjunction is the conjunction of the complements.It is used to calculate the "opposite" of a solution.
In optimization algorithms, a point and its opposite define an axis that bisects the solution space.The point lies on one side of the space, and its opposite is equidistant from the centroid but on the opposite side of the space.This geometric interpretation offers the key insight of OBL: for each point considered by the algorithm, there exists an unexplored point that is as far from the centroid as the original point but in the opposite direction.If a solution is a point in the search space, its opposite is another point in the space such that the line connecting the solution and its opposite passes through the centroid of the search space.The primary advantage of considering both a solution and its opposite is that it can double the useful information obtained from each fitness evaluation.Furthermore, because the opposite points are spread throughout the solution space, OBL can widen the search scope, enhancing exploration capabilities and accelerating convergence to the global optimum.Therefore, the proposed algorithm can explore the search space more effectively by evaluating both a solution and its opposite.The OBL is a strategy used to enhance the performance of optimization algorithms by exploring the search space more efficiently.The central concept of OBL is the simultaneous consideration of an estimate (a solution) and its opposition during the search process.In mathematical terms, for a problem defined in the dim-dimensional real space, a point X w in the search space is represented as a vector of dim real numbers.Given a real-valued solution " X w " in the interval [lb, ub] , the opposite solution " X w,opp " is calculated as follows.
where lb and ub represent the lower and upper bounds of the search space, respectively.The fitness of both X w and X w,opp are evaluated, and the point with the better fitness is selected, i.e., iff (X w,opp ) < f (X w ), thenX w = X w,opp .
The main idea of this strategy is to consider the current estimate and its opposite to obtain a better approximation for the global optimal solution.In other words, for each candidate solution, the algorithm generates its opposite solution and evaluates both.The algorithm tries to keep the one with better fitness for the next generation.This can be especially useful in the initial stages of the algorithm's run, where it can significantly increase the convergence speed.The reason is that, with OBL, the algorithm can simultaneously consider two "opposite" points in the solution space.This is akin to exploring two different directions simultaneously, potentially leading to a broader and more efficient search.This OBL strategy is incorporated into EDO during the initialization and mutation phases.As discussed, during initialization, each randomly generated solution X w is accompanied by its opposite X w,opp , and the one with better fitness is selected for the initial population.An offspring solution is created, and its opposite is generated during mutation.Again, the one with the better fitness is selected to replace its parent.This way, in each iteration, the algorithm explores the vicinity of the solutions and their opposite points in the search space.This strategy significantly enhances the exploratory capability of the algorithm, making it more robust and potentially faster in finding the global optimal solution.By applying this strategy, the OBEDO can benefit from an enhanced exploration capability, potentially improving its effectiveness and efficiency in (28) Initialize the parameters, such as , , bounds, etc. Initialize the random population and generate the random population solutions using Eq.14.
Generate the opposite population solutions using Eq.29.

End if
Construct the memory matrix.While ( < ) do Define the matrix of size ∈ Compute using Eq. 15 Define all the adaptive parameters , , , , Update the position using Eq.16 (Exploitation phase) else Update the position using Eq.23 (Exploration phase)

End if
Copy the matrix to the memoryless matrix.Update the solution comparing the current fitness and updated fitness.
End while Return: BestPosition and the BestFitness.

Algorithm 2. Pseudocode of the proposed OBEDO algorithm.
The OBEDO algorithm, as discussed, likely integrates two key concepts: OBL and EDO algorithm.OBL is a concept used in optimization algorithms to speed up the convergence rate.It works on the principle that considering a candidate solution and its opposite can provide a better approximation of the global optimum.In OBL, for every estimated solution, the opposite solution is also considered.This means if a solution is at a point x in the search space, its opposite −x is also evaluated.This approach increases the chances of finding better solutions in fewer iterations as it explores the search space more effectively.When integrated into the OBEDO framework, OBL could help in quickly identifying more promising regions of the search space for parameter estimation, thereby improving the efficiency and accuracy of the optimization process.OBL enhances the global search capability, ensuring diverse and comprehensive exploration, while the exponential distribution method fine-tunes the search, allowing for efficient exploitation of promising areas.This synergy could make OBEDO particularly effective in navigating complex, high-dimensional search spaces typical in solar PV model parameter estimation, where traditional methods might struggle due to local optima or slow convergence.The effectiveness of OBEDO in outperforming other methods could stem from its ability to balance exploration and exploitation in the optimization process.This balance is crucial in parameter estimation problems, where finding the global optimum in a complex search space is essential for accurate and reliable results.

Complexity of the proposed OBEDO algorithm
The algorithm creates an initial population of solutions X w i .The complexity of the initialization phase should be O(N p × dim) , where N p is the number of individuals, and dim is the dimension of the solution space.The algorithm performs a series of computations for each individual in the population ( N p ), and this is done for each iteration.Thus, the complexity within the loop could be considered as O(N p × dim) .This includes fitness computation, memory updates, and solution modifications.Assuming the complexity of the fitness function is O(dim) , this would still give a complexity of O(N p × dim) for the whole iteration.The sorting operations would add a complexity of O(N p × log(N p )) per iteration.However, since N p × log(N p ) is smaller than N p × dim for large enough N p and dim , ignore it in the final complexity assessment.Therefore, the time complexity of the OBEDO can be approximated as O(N p × dim × T) , which is the same as the original EDO algorithm.

Results and discussions
In this section of the paper, we present the experimental results obtained to assess the performance of the proposed OBEDO algorithm.The OBEDO determines the parameters of four benchmarked PV cell/module models: RTC France silicon cell, PVM752 GaAs PV cell, Photowatt PWM201 module, commercial Sharp ND-R250A5 PV module, and commercial SM55 PV module.To evaluate the effectiveness of the OBEDO, it is applied to three PV models, namely the SDM, DDM, and TDM.For each model, the OBEDO and other seven algorithms, such as opposition-based GBO (OBGBO), opposition-based marine predator algorithm (OBMPA), MGTO, ADHHO, IAOA, HDE, and original EDO, are employed to obtain the optimal parameters that best fit the experimental data.The application of the proposed OBEDO and other algorithms to estimate PV parameters is pictorially represented in Fig. 5.
For the RTC France silicon cell with a diameter of 57 mm, experimental data is collected under specific conditions, with a solar irradiance of 1000 W/m 2 and a temperature of 33 °C.A total of 26 sets of current and voltage measurements are gathered to determine the cell variables by optimizing the RMSE, a common metric for evaluating the accuracy of models.The goal is to find the best-fitting parameters that minimize the RMSE and, thus, enhance the accuracy of the PV cell model.For the PVM752 GaAs PV cell, the experimental data is gathered under 1000W/m 2 irradiation and 25 °C temperature.Similarly, for the Photowatt PWM201 module, experimental data is gathered under specific conditions with a solar insolation of 1000 W/m 2 and a temperature of 45 °C.For the commercial Sharp ND-R250A5 PV module, the experimental data is gathered under 1040W/ m 2 irradiation and 59 °C temperature.In addition, the data for the commercial SM55 PV module under different operating conditions are also collected.To ensure realistic parameter estimates, lb and ub limits for the cells and the modules parameters are provided in Table 2.These limits constrain the optimization process during the execution of all algorithms and ensure that the resulting parameters fall within physically meaningful ranges.
The simulation tests were conducted on a laptop running Windows 11 with specific hardware specifications.The laptop has an Intel(R) Core (TM) i5-10300H CPU operating at a clock frequency of 2.44 GHz and 8 GB of RAM.For the simulation tests, this study considered a population size of 40, and the maximum number of iterations was set to 1000 for all the PV models under consideration.Additionally, other control parameters for all the algorithms utilized in the experiments are documented in Table 3, providing transparency and reproducibility of the study.To ensure a fair comparison among all the selected algorithms, each method was executed 30 times independently.This repetition helps account for any potential variability in the results and allows us to draw robust conclusions regarding the algorithms' performance.This study employed MATLAB R2020b software to conduct the simulation tests.By running the experiments on the specified laptop and using the common platform of MATLAB R2020b, this study ensures consistency and comparability of the results across all algorithms and  PV models.This systematic approach allows us to make well-informed evaluations of the algorithms' efficiency and effectiveness in optimizing the PV model parameters.
In this study, the performance comparison among all the selected algorithms is carried out using various statistical measures and performance metrics to assess their effectiveness in optimizing the PV model parameters.The statistical measures considered include Minimum (Min), Maximum (Max), Mean, and Standard Deviation (STD).These measures help us understand the range and distribution of the results obtained by each algorithm.On the other hand, performance metrics provide specific quantitative evaluations of the algorithms' accuracy in predicting the PV model parameters.The metrics used in this comparison are Relative Error (RE), Integral Absolute Error (IAE), and Root Mean Square Error (RMSE).The runtime (RT) of each algorithm is also analyzed, which represents the average running time of all 30 independent runs for each case study.This information is crucial for understanding the computational efficiency of the algorithms and their feasibility for practical applications.The metrics RE and IAE are computed as follows: (i) IAE represents the integral absolute error for a particular trial, which is calculated as the absolute difference between the estimated current I es and the measured current I ex value for that trial, and (ii) RE denotes the relative error, which quantifies the percent- age difference between the estimated current and the measured current value.It is computed as the ratio of the difference between the measured and estimated current to the measured current value.The IAE is designed   www.nature.com/scientificreports/statistical measures, performance metrics, runtime analysis, and Friedman's Ranking Test, the study comprehensively evaluates the selected algorithms' performance and provides valuable insights into their capabilities in accurately estimating the PV model parameters.Such detailed comparisons are essential for researchers and practitioners to decide on the most suitable algorithm for specific PV modelling applications.This sub-section details the results obtained by the proposed algorithm and other algorithms for scenario 1, i.e., SDM, DDM, and TDM of the RTC France Si PV cell.The bounds for all the PV models are provided in Table 2.
As deliberated earlier, all algorithms are executed 30 times for a fair comparison.Tables 4, 5, and 6 display a comprehensive overview of the five, seven, and nine parameters estimated through the SDM, DDM, and TDM employment pertaining to the RTC of France Si solar cell.It is imperative to emphasize that these parameters were determined within the precise confines, as elucidated in Table 2.A graphical representation of the I-V (Current-Voltage) and P-V (Power-Voltage) characteristic curves pertinent to the SDM, DDM, and TDM of the RTC France Si solar cell obtained by all algorithms have been expounded in Figs. 6, 7, and 8.These elucidate the outcomes engendered by the proposed OBEDO that was instituted as part of this study.It is crucial to underline that the quantitative values encapsulating the IAE and the RE have been meticulously documented in Table 7.To visually portray the paramount significance of the RE and IAE values ascertained through the instrumentality of the OBEDO and other algorithms, this study has deviously created error graphs as showcased in Figs. 9, 10, and 11.A cursory examination of Figs. 6, 7, 8, 9, 10, and 11 unfolds a noticeably harmonious alignment between the projected estimated and empirically derived experimental curves.This alignment strengthens the inference that the curve fitting has yielded a commendable equivalence.To offer a comprehensive assessment, meticulous cataloguing of statistical indicators including, but not limited to, Min, Mean, Max, RT, and STD values has been meticulously undertaken.These values have been scrupulously logged within Table 8 to facilitate in-depth data comprehension.Furthermore, it is noteworthy that salient achievements, exemplified by the most remarkable outcomes gleaned from the diverse tables, have been judiciously highlighted through the distinctive formatting of boldfaces.
The IAE value resulting from using the OBEDO in the context of the SDM, DDM, and TDM for scenario 1 is notably below the threshold of 3.01E−03.This indicates a remarkably precise alignment between the predicted and observed values.Similarly, the RE value, which gauges the degree of dissimilarity between the predicted and actual values, demonstrates significant conformity as it remains under the stringent limit of 1.51E−01.A broader assessment of the performance metrics sheds light on the collective behaviour for scenario 1 under the purview of the OBEDO algorithm.The average RE and IAE values of SDM, DDM, and TDM are 4.74E−03, 4.50E−03, 4.51E−03, 8.28E−04, 8.19E−04, and 8.18E−04respectively, further reinforcing the efficacy of the optimization algorithm in achieving a robust convergence between the simulated and empirical outcomes.The comprehensive data showcased in Table 5 and Figs. 9, 10, and 11 corroborates the capability of the OBEDO to precisely extract the intricate characteristics of the PV cell within the specific framework.This affirmation of accuracy underscores the algorithm's competence in capturing the nuanced behaviour and performance attributes of the PV cell, thereby contributing to a deeper and more comprehensive understanding of its indefinite characteristics.Table 8 summarizes the performance of various algorithms across three models: SDM, DDM, and TDM.The metrics assessed include Min, Max, Mean, Median, STD, and RT for each algorithm.Across all three models, the OBEDO algorithm consistently achieves low RMSE values, indicative of accurate predictions.Notably, it maintains a narrow RMSE range, resulting in stable and reliable outcomes.The EDO algorithm displays more variability in RMSE, particularly in the DDM and TDM models, where it reaches higher maximum RMSE values.This suggests that EDO performs well on average but may produce less precise predictions occasionally.ADHHO demonstrates generally consistent performance, maintaining moderate RMSE values and standard deviations.Its RMSE range remains relatively stable across the different models.OBMPA consistently produces moderate RMSE values with low standard deviations, indicating dependable performance.Its RMSE range is also relatively stable across models.MGTO consistently achieves low RMSE values with minimal variability, suggesting accurate and reliable predictions.Notably, its runtime is consistent across models.IAOA displays wider variability in RMSE, particularly in the SDM and DDM models, where it reaches higher maximum values.While its mean RMSE is reasonable, the variability suggests some potential for less accurate predictions.HDE generally maintains low and consistent RMSE values, with minor fluctuations.Its RMSE range remains stable across models, indicating dependable performance.OBGBO stands out with extremely low RMSE values, particularly in the SDM model.However, it comes with higher runtimes, especially in the TDM model.The OBEDO, OBGBO, and MGTO algorithms consistently exhibit accurate and reliable performance across the three models.OBMPA, EDO, ADHHO, IAOA, and HDE also show reasonable performance but with some variability.OBGBO demonstrates good accuracy, though at the cost of longer runtimes.The progression trajectories, characterizing the convergence behaviour, intricately unravel the dynamics of the SDM, DDM, and TDM for Scenario 1.These complex patterns of convergence are portrayed in the comprehensive Fig. 12.As the gaze fixes upon this graphical representation, the remarkably accelerated convergence of the OBEDO becomes noticeably apparent.This swift trajectory toward convergence, markedly outpacing its algorithmic counterparts, is a testament to the superior efficiency of the OBEDO.The unfolding narrative takes an artistic turn as we venture into Fig. 13, which, rather than opting for the conventional boxplot approach, ingeniously employs the visually captivating violin plot technique.Within this illustration, the distinct strategies outlined for the SDM, DDM, and TDM emerge as leading roles, each contributing its unique melodic note to the symphony of data visualization.In the symphony of violin plots, a striking revelation comes to the fore-OBEDO's reliability emerges as a true virtuoso.A harmonious refrain of low STD values resonates across all the photovoltaic models in Scenario 1.
This explicit stability display strengthens the OBEDO's position as an unmatched performer, standing high amidst its algorithmic companions.Interestingly, the plot's narrative takes a dramatic twist as a few algorithms stumble upon the snares of local minima, their quest for optimal outcomes thwarted by the intricate labyrinth of possibilities.Yet, amidst this unfolding drama emerges the OBEDO as the triumphant protagonist.It directs the complex optimization landscape, ascending to the peak of achievement with incomparable accuracy.

Scenario 2-PVM752 GaAs PV Cell
This sub-section details the results obtained by the proposed algorithm and other algorithms for scenario 2, i.e., SDM, DDM, and TDM of the PVM752 GaAs PV cell.The bounds for all the PV models are provided in Table 2.   13 to facilitate in-depth data comprehension.
The IAE value resulting from using the OBEDO in the context of the SDM, DDM, and TDM for scenario 2 is notably below the threshold of 9.01E−04.This indicates a remarkably precise alignment between the predicted and observed values.Similarly, the RE value, which gauges the degree of dissimilarity between the predicted and actual values, demonstrates significant conformity as it remains under the stringent limit of 1.51E−02.The average RE and IAE values of SDM, DDM, and TDM are − 7.46E−03, − 5.84E−03, − 6.67E−03, 2.02E−04, 1.63E−04, and 1.36E−04, respectively, further reinforcing the efficacy of the optimization algorithm in achieving a robust convergence between the simulated and empirical outcomes.The comprehensive data showcased in Table 12 and Figs.17, 18, and 19 corroborates the capability of the OBEDO to precisely extract the intricate characteristics of the PV cell within the specific framework.This affirmation of accuracy underscores the algorithm's competence in capturing the nuanced behaviour and performance attributes of the PV cell, thereby contributing to a deeper and more comprehensive understanding of its indefinite characteristics.
Table 13 summarizes the performance of various algorithms across three models: SDM, DDM, and TDM.Across all the PV models, the algorithms exhibit notable variations in their capabilities.The OBEDO algorithm consistently demonstrates commendable accuracy, as reflected by its frequently low IAE values, ensuring precise predictions.Contrasting this, the EDO algorithm presents a broader spectrum of RMSE values, potentially indicating instances where predictions may be less precise.The ADHHO algorithm maintains consistent and moderate RMSE values, reflecting stable performance across the PV models.Similarly, the OBMPA algorithm displays uniform RMSE values, suggesting reliable outcomes across different scenarios.The MGTO algorithm consistently produces low RMSE values with minimal variability, indicating its potential for reliable predictions.Meanwhile, the IAOA algorithm exhibits more variability in RMSE, especially in the SDM and DDM models.Remarkably, the HDE algorithm consistently maintains low and consistent IAE values across all PV models.The OBGBO algorithm stands out for its exceptional accuracy, albeit with longer runtimes.In summary, the algorithms' performance reveals a diverse landscape of accuracy, convergence, and efficiency.While some algorithms shine in specific models, the OBEDO, MGTO, and HDE algorithms consistently exhibit robust and dependable performance across various photovoltaic models.The findings highlight the crucial interplay between accuracy and efficiency, underscoring the necessity for a comprehensive evaluation of optimization algorithms in photovoltaic modelling.The IAOA and OBMPA failed to trace the I-V and curves for the SDM because both get trapped by the local optima.
The convergence dynamics of the SDM, DDM, and TDM for Scenario 2 are meticulously unveiled through intricate progression trajectories, portrayed vividly in Fig. 20.Upon inspecting this visual depiction, the rapid and impressive convergence of OBEDO becomes clearly evident.This accelerated trajectory surpasses its algorithmic counterparts, underscoring OBEDO's exceptional efficiency.The narrative takes an artistic turn with Fig. 21, which employs the visually captivating violin plot.The strategies defined for SDM, DDM, and TDM play leading roles, each contributing a distinct note to the symphony of data visualization.Within the realm of the symphony of violin plots, a remarkable revelation unfolds-OBEDO emerges as a virtuoso of reliability.The attentive observer discerns the delicate undercurrents of STD values.The low STD values resonate across all photovoltaic models in Scenario 2. This steadfast stability further solidifies OBEDO's stature as an unmatched performer, towering among its algorithmic peers.convergence between the simulated and empirical outcomes.The comprehensive data is showcased in Table 17 and Figs. 25, 26, and 27 corroborates the capability of the OBEDO to precisely extract the intricate characteristics of the PV module within the specific framework.This affirmation of accuracy underscores the algorithm's competence in capturing the nuanced behaviour and performance attributes of the PV module, thereby contributing to a deeper and more comprehensive understanding of its indefinite characteristics.
Table 18 provides a comprehensive overview of algorithmic performance across three photovoltaic models: SDM, DDM, and TDM.Key metrics are presented for each algorithm within these models.Across the PV models, a diverse array of algorithmic behaviour is evident.OBEDO consistently demonstrates remarkable accuracy with low RMSE values, indicating precise predictions.EDO shows broader variability in RMSE values,  MGTO and IAOA present varying precision, while HDE demonstrates consistent accuracy.OBGBO maintains its performance pattern with higher runtimes.Overall, OBEDO, MGTO, and HDE consistently exhibit robust performance across various PV models, with OBEDO excelling in accuracy and efficiency.The findings emphasize the intricate interplay between accuracy and runtime, underscoring the significance of a comprehensive evaluation of algorithmic efficacy in PV modelling.The convergence dynamics of the SDM, DDM, and TDM for Scenario 3 are meticulously unveiled through intricate progression trajectories, portrayed vividly in Fig. 28.Upon inspecting this visual depiction, the rapid and impressive convergence of OBEDO becomes clearly evident.This accelerated trajectory surpasses its algorithmic counterparts, underscoring OBEDO's exceptional efficiency.The narrative takes an artistic turn with Fig. 29, which employs the visually captivating violin plot.
The strategies defined for SDM, DDM, and TDM play leading roles, each contributing a distinct note to the symphony of data visualization.Within the realm of the symphony of violin plots, a remarkable revelation unfolds-OBEDO emerges as a virtuoso of reliability.The attentive observer discerns the delicate undercurrents of STD values.The low STD values resonate across all photovoltaic models in Scenario 3.This steadfast stability further solidifies OBEDO's stature as an unmatched performer, towering among its algorithmic peers.

Scenario 4: commercial sharp ND− R250A5 polycrystalline PV module
The I-V characteristic of this panel has been measured at a temperature of 59 °C and an irradiation of 1040 W/m 2 .The results obtained for SDM of the commercial Sharp ND-R250A5 polycrystalline PV module are discussed in this sub-section.This PV panel consists of 60 polycrystalline silicon PV cells that are connected in series.Table 19 presents the optimized parameters of the commercial Sharp ND-R250A5 polycrystalline PV module by all selected algorithms.Table 20 presents the recorded data points, estimated data points using the OBEDO, and the corresponding errors.The average values for IAE and RE are stated as 9.41E−03 and − 3.93e−04.Table 2 displays the boundary values, and Table 19 displays the estimated optimal parameters of the PV module calculated by the OBEDO.The statistical measures for the SDM of this PV panel are presented in Table 21.The optimal RMSE is 1.12E−02, while the IAE is 9.41e−03.The results of five parameters derived by OBEDO were given in Table 19 and were compared with the other statE−of-thE−art algorithms.When comparing the outputs of the SDM, it is evident that the SDM obtained by OBEDO has lower RMSE values compared to the other algorithms.The error characteristics of this case are visualized in Fig. 30.Figure 30a illustrates the IAE with respect to the measured voltage, and Fig. 30b illustrates the RE with respect to the measured voltage.The measured and estimated I-V characteristics for the SDM of this PV panel are depicted in Fig. 31a.The P-V characteristics of this PV panel are shown in Fig. 31b, both measured and computed.The convergence of the objective function and violin plot analysis for OBEDO and seven other optimisation strategies is shown in Fig. 32.The convergence speed of the OBEDO surpasses that of the other algorithms, as demonstrated by Fig. 32a.The reliability of all the algorithms is visualized in Fig. 32b.

Scenario 5-commercial SM55 PV module
This sub-section details the results obtained by the proposed algorithm and other algorithms for scenario 5, i.e., the SDM of the SM55 PV module.The bounds for this scenario are provided in Table 2. To extract the parameters of scenario 5, a comprehensive analysis is conducted using experimental I-V samples.These samples are gathered under specific conditions, encompassing constant irradiance (G) of 1000 W/m 2 and varying temperatures (T)   In the context of this study, the variables T and G hold significance as they denote the prevailing temperature and irradiation levels.The parameter α signifies the temperature coefficient, while T STC represents the temperature at STC, and G STC corresponds to the irradiance at STC.Additionally, I (sc(STC)) represents the short-circuit current (30) Table 16.Parameters estimated for Scenario 3 (TDM).at STC.Through this meticulous examination, insights are gained into the behaviour and adaptability of the PV module model, facilitating a nuanced understanding of its performance under varying temperature and irradiance conditions.This analysis contributes to a comprehensive comprehension of the model's capabilities, aiding in optimizing and deploying photovoltaic systems across a spectrum of real-world scenarios.Within this study, Tables 22 and 23 play a pivotal role in offering a comprehensive insight into the parameters estimated and optimized using the OBEDO, alongside several alternative methodologies employed to investigate the PV module model.Delving into the specifics, Table 22 meticulously catalogues the estimated variables associated with the SM55 photovoltaic module.The entries in this table pertain to scenarios where the irradiance remains constant at 1000 W/m 2 while the temperature fluctuates across three distinct levels: 25 °C, 50 °C, and 75 °C.Meanwhile, Table 23 presents a similar compilation of estimated parameters for the SM55 PV module under differing conditions.In this instance, the temperature is maintained at a constant 25 °C while the irradiance levels are varied, spanning the range of 1000-200 W/m 2 .
The outcomes in Tables 22 and 23 witness the superior predictive prowess of the proposed OBEDO algorithm.These assertions are substantiated by observing significant RMSE values.This quantifiable measure of accuracy underscores the algorithm's efficacy in determining and predicting the parameters intrinsic to a PV system with enhanced precision.Augmenting the tabular data, Figs.33 and 34 provide a graphical representation of the I-V and P-V curves, depicting the behaviour of the SM55 PV module as extracted through the OBEDO and other algorithms.These curves offer a tangible means to observe the alignment and coherence between the curves generated from the algorithm's predicted data and the actual empirical data.This visual correlation elucidates  how much the algorithmic predictions harmonize with real-world observations, enhancing our understanding of the model's performance and predictive capabilities.
In order to conduct a more comprehensive evaluation of the performance of the newly introduced OBEDO, a detailed analysis of its statistical measures was carried out.This analysis specifically focused on Scenario 5 and involved the computation of various key metrics, including the Min, Max, Mean, Median and STD values.These statistical results were subsequently organized and presented in a tabular format, with Tables 24 and 25 dedicated to showcasing these values.The investigation was undertaken across various environmental factors, encompassing different temperature settings and varying irradiation conditions.The purpose of this exploration was to measure the algorithm's robustness and effectiveness under a variety of scenarios.Notably, the outcomes  yielded by the OBEDO exhibited remarkable performance in terms of the computed statistical measures, specifically the Min, Max, Mean, and STD values.This notable achievement can be attributed to incorporating a unique enhancement technique, namely integrating OBL within the algorithmic framework.By infusing the algorithm with the principles of OBL, the inherent capacity of the OBEDO to explore and exploit the solution space was significantly augmented.This enhancement played a crucial role in the algorithm's ability to search for optimal PV module parameters under varying conditions efficiently.The synergistic effect of the OBL strategy, alongside other optimization techniques embedded within the OBEDO, contributed to the algorithm's exceptional performance across different environmental contexts.In summary, the meticulous evaluation of the proposed OBEDO's performance in Scenario 5, through the assessment of statistical measures across distinct temperatures and irradiation conditions, underscored its effectiveness in achieving outstanding results.This success can be attributed to the successful integration of OBL, which bolstered both exploration and exploitation capabilities, ultimately identifying optimal PV module parameters across various conditions.However, it is worth noting that the RT values acquired through utilizing the proposed algorithm across all operational conditions exhibit a marginal increase compared to those attained by the fundamental EDO.This observed difference in RT values can be attributed to amalgamating two distinct optimization strategies within the proposed algorithmic framework.Despite the slightly elevated RT values, the overall performance of the proposed algorithm surpasses that of the basic EDO by a notable margin, showcasing an improvement of over 75%.This significant enhancement can be attributed to the synergistic effect generated by the fusion of these diverse optimization strategies.The results obtained from the evaluation demonstrate a remarkable alignment between the I-V and P-V characteristics deduced from the estimated parameters and the experimental data, even    under varying temperature and irradiance conditions.This alignment underscores the precision of the proposed algorithm in capturing the intricate nuances of the photovoltaic system's behaviour.Furthermore, the outcomes of the experimentation reveal an additional advantage of the OBEDO, namely its ability to attain a lower RMSE value.This improvement in accuracy highlights the algorithm's efficacy in modelling and predicting the behaviour of the photovoltaic module, thereby facilitating more reliable parameter estimations.The comprehensive discussions and analyses led to the clear conclusion that the efficiency of the OBEDO remains robust and dependable when confronted with dynamic shifts in environmental conditions.The amalgamation of optimization strategies, while contributing to a slight increase in RT values, offers substantial gains in performance, as evident from the     substantial enhancement over the basic EDO.The alignment between simulated I-V and P-V characteristics and experimental data, coupled with the lower RMSE value achieved, attests to the algorithm's capability to handle the details of dynamic environmental variations.

Statistical performance
Milton Friedman, a renowned figure in statistics, is credited with conceiving a significant non-parametric statistical test known as the Friedman ranking test (FRT).This test holds a distinct purpose in the realm of statistical analysis-it serves as a valuable tool for detecting variations in treatments across multiple experimental runs, akin to the parameterized repeated measures ANOVA commonly employed in research.A comprehensive application of the FRT was conducted to establish the superiority of the OBEDO in terms of overall performance.In this endeavour, the objective was to substantiate that the OBEDO exhibits superior efficacy in generating aggregate results compared to a range of other algorithms.Among the algorithms subjected to analysis were not only the OBEDO and EDO but also an assortment of others, namely the ADHHO, OBMPA, MGTO, IAOA, HDE, and OBGBO.The process involved meticulously examining the FRT outcomes, which were carefully obtained by considering each algorithm's RMSE standard deviation values.These values, extracted and compiled in Table 26, provided a foundation for an extensive comparative analysis of the various algorithms under investigation.The resultant observations from this comparison laid bare a noteworthy revelation -the OBEDO emerged as a standout performer.By assessing the mean FRT values across four distinct case studies and comparing them against the benchmark set by preceding algorithms, the authors concluded that the recommended OBEDO showcases an exceptional ability to manage diverse scenarios with the lowest mean FRT values.In light of the compelling evidence amassed through this comprehensive analysis, it becomes evident that the OBEDO outshines its counterparts.The OBEDO's consistent superiority across multiple case studies firmly establishes its prowess as a premier algorithmic solution, further reinforcing its position as a robust and proficient tool for addressing various optimization challenges.Table 26 presents a comparative analysis using the FRT across various scenarios and optimization models.The FRT aims to determine the relative performance of different optimization algorithms.The models evaluated Notably, OBEDO maintains the top rank with the lowest mean FRT value, emphasizing its consistent superiority.In summary, the FRT analysis across multiple scenarios and models underscores the exceptional performance of the OBEDO algorithm, which consistently achieves the lowest FRT values and the highest ranking.This suggests that OBEDO outperforms the other optimization models across diverse conditions and demonstrates its robustness and efficacy in optimization tasks.

Conclusions and future extension
The proposed OBEDO has demonstrated its prowess in efficiently, accurately, and rapidly extracting undefined parameters from diverse PV system models, as evidenced by the proposed study.The comprehensive analysis undertaken further substantiates the efficacy and potential of this framework in extrapolating parameters for commercial PV models.The evaluation sought to ascertain if the proposed methodology could successfully deduce the parameters of such intricate models, and the obtained results provide compelling evidence of the OBEDO's capability to precisely and effectively determine optimal parameters.Furthermore, it is noteworthy that the TDM exhibits the lowest model features in terms of RMSE, IAE, and RE, underscoring its accuracy.Additionally, the convergence speed of the OBEDO outperforms its counterparts across all test cases, indicating its superior efficiency.These attributes collectively position the proposed OBEDO as a cutting-edge, competitive, and advanced approach for parameter identification in PV models.The OBEDO attains stability and diversity, achieved through meticulous parameter tuning, which positions it favourably compared to its basic version.Nonetheless, as per the NFL theorem, it is crucial to recognize that while the OBEDO exhibits enhanced search results for parameter identification problems, its success might vary in problems with distinct characteristics.As a result, ongoing research and development are necessary to tailor meta-heuristic search approaches to specific problem domains.
In forthcoming research, the proposed algorithm can be evaluated within contemporary techniques for PV parameter identification and subsequently extended to real-world scenarios for various PV cells/modules.Moreover, the versatile applicability of the OBEDO can be harnessed to address challenging problems spanning water pollution prediction, disease diagnosis, binary optimization, feature selection/extraction, image processing tasks, scheduling optimization, wireless sensor networks, and medical image categorization.This underscores the algorithm's potential as a multifaceted solution with far-reaching implications across diverse domains.

Figure 1 .
Figure 1.SDM of the photovoltaic cell.

Figure 2 .
Figure 2. DDM of the photovoltaic cell.

Figure 3 .
Figure 3. TDM of the photovoltaic cell.

Figure 5 .
Figure 5. Application of the OBEDO and other algorithms for parameter estimation.

Scenario 3 -
www.nature.com/scientificreports/Photowatt PWP-201 PV module This sub-section details the results obtained by the proposed algorithm and other algorithms for scenario 3, i.e., SDM, DDM, and TDM of the Photowatt PWP-201 PV module.The bounds for all the PV models are provided in Table 2. Tables 14, 15, and 16 display a comprehensive overview of the five, seven, and nine parameters estimated through the SDM, DDM, and TDM employment pertaining to the Photowatt PWP-201 PV module.A graphical representation of the I−V and P-V characteristic curves pertinent to the SDM, DDM, and TDM of the Photowatt PWP-201 PV module obtained by all algorithms has been expounded in Figs.22, 23, and 24.The quantitative values summarizing the IAE and the RE must be meticulously documented in Table 17.To visually portray the paramount significance of the RE and IAE values ascertained through the instrumentality of the OBEDO and other algorithms, this study has deviously created error graphs as showcased in Figs. 25, 26, and 27.A cursory examination of Figs.22, 23, 24, 25, 26, and 27 unfolds a noticeably perfect alignment between the estimated and experimental curves.This alignment strengthens the inference that the curve fitting has yielded a worthy equivalence.To offer a comprehensive assessment, meticulous cataloguing of statistical indicators has been meticulously undertaken.These values have been carefully logged within Table18to facilitate in-depth data comprehension.The IAE value resulting from using the OBEDO in the context of the SDM, DDM, and TDM for scenario 2 is notably below the threshold of 5.12E−04.This indicates a remarkably precise alignment between the predicted and observed values.Similarly, the RE value, which gauges the degree of dissimilarity between the predicted and actual values, demonstrates significant conformity as it remains under the stringent limit of 2.41E−02.The average RE and IAE values of SDM, DDM, and TDM are − 3.253E−04, 3.253E−04, 3.253E−04, 1.957E−03, 1.957E−03, and 1.957E−03, respectively, further reinforcing the efficacy of the optimization algorithm in achieving a robust
of 25 °C, 50 °C, and 75 °C, as well as a constant temperature of 25 °C and diverse irradiation levels including 1000 W/m 2 , 800 W/m 2 , 600 W/m 2 , 400 W/m 2 , and 200 W/m 2 .The evaluation process encompasses various statistical measures and other pertinent performance indicators.Equation (30) is pivotal in this context, as it calculates the specific photovoltaic modules short-circuit current I sc across various operational scenarios.

Figure 33 .
Figure 33.Characteristics of the SM55 PV module under different temperature conditions obtained by all algorithms; (a) I-V curves, (b) P-V curves.

Figure 34 .
Figure 34.Characteristics of the SM55 PV module under different irradiance conditions obtained by all algorithms; (a) I-V curves, (b) P-V curves.
2)I pv = I ph − I sd exp V pv + I pv R se nV t − 1 − V pv + I pv R se R sh 2, . . ., N p andj = 1, 2, . . ., dim solving complex optimization problems.The pseudocode of the proposed OBEDO algorithm is presented in Algorithm 2. The flowchart of the OBEDO is shown in Fig.4.

Table 3 .
Parameters of all the algorithms.

Table 5 .
Parameters estimated for Scenario 1 (DDM).Significant values are in [bold].Vol.:(0123456789)Scientific Reports | (2024) 14:528 | https://doi.org/10.1038/s41598-023-50890-ywww.nature.com/scientificreports/ to penalize errors equally, regardless of the direction, while the RE provides insight into the magnitude of the absolute error with respect to the measured data.When an exact measurement is unavailable, using the measured value to calculate relative inaccuracy is common.Finally, a Friedman's Ranking Test (FRT) is performed to validate the statistical significance of the results.This test allows for comparing multiple algorithms across different metrics and identifies whether any algorithm significantly outperforms the others.By employing these

Table 7 .
RE and IAE values attained by OBEDO for Scenario 1.

Table 8 .
Performance comparison for Scenario 1. Significant values are in [bold].

Table 12 .
RE and IAE values attained by OBEDO for Scenario 2.

Table 13 .
Performance Comparison for Scenario 2. Significant values are in [bold].

Table 17 .
RE and IAE values attained by OBEDO for Scenario 3.

Table 18 .
Performance Comparison for Scenario 3. Significant values are in [bold].

Table 19 .
Parameters estimated for Scenario 4. Significant values are in [bold].

Table 20 .
RE and IAE values attained by OBEDO for Scenario 4.

Table 21 .
Performance comparison for Scenario 4. Significant values are in [bold].

Table 22 .
Parameters estimated by the algorithms (Scenario 4) under various temperature conditions.Significant values are in[bold].SDM, DDM, TDM), the FRT values indicate the algorithms' performance.The OBEDO algorithm achieved the lowest FRT value, followed by MGTO and OBGBO.Similar trends are observed in Scenario 2, which considers another set of data models, with OBEDO outperforming other models.Scenario 3 explores algorithms' performance across diverse conditions.Once again, OBEDO takes the lead in achieving the lowest FRT value, followed by ADHHO and IAOA.Scenario 4 explores algorithms' performance across diverse conditions.Once again, OBEDO takes the lead in achieving the lowest FRT value, followed by OBGBO and MGTO.In the fifth scenario, different temperature and irradiance conditions are analyzed.OBEDO consistently demonstrates the most favourable performance across varying environmental conditions.The mean FRT provides the average FRT values for each model, and the rank presents their ranking based on performance.

Table 23 .
Parameters estimated by the algorithms (Scenario 4) under irradiation conditions.Significant values are in [bold].

Table 24 .
Statistical analysis of Scenario 4 at different temperatures.Significant values are in [bold].

Table 25 .
Statistical analysis of Scenario 4 at different irradiances.Significant values are in [bold].

Table 26 .
Statistical performance comparison.Significant values are in [bold].